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(54, Automated real-time analysis of nucleic acid amplification 

(57) A method is described for analyzing a sample 
for the presence of a nucleic acid wherein the sample is 
amplified, preferably using PCR, in the presence of a 
fluorescent probe capable of detecting the presence of 
the nucleic acid sample. A baseline region is deter- 
mined by comparing the fluorescence at various ampli- 
fication cycles, and the fluorescence at each of various 
amplification cycles is compared to the baseline region 
to determine whether the fluorescence measurements 
fall outside of that baseline region. 
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Description 

Beidoj the Invention 

[0001] The present invention relates to a method of 
analyzing a sample for the presence of a nucleic acid. 
More particularly, the present invention is directed to an 
automated method for detecting and reporting the pres- 
ence of a predetermined nucleic acid in a sample using 
polymerase chain reaction and a fluorescent detecting 
entity. 

Rark gmiind and Summary of the Invention 

[0002] Amplification of DNA by polymerase chain 
reaction (PCR) is a technique fundamental to molecular 
biology. Nucleic acid analysis by PCR requires sample 
preparation, amplification, and product analysis. 
Although these steps are usually performed sequen- 
tially, amplification and analysis can occur simultane- 
ously. DNA dyes or fluorescent probes can be added to 
the PCR mixture before amplification and used to ana- 
lyze PCR products during amplification. Sample analy- 
sis occurs concurrently with amplification in the same 
tube within the same instrument. This, combined 
approach decreases sample handling, saves time, and 
greatly reduces the risk of product contamination for 
subsequent reactions, as there is no need to remove the 
samples from their closed containers for further analy- 
sis. The concept of combining amplification with product 
analysis has become known as 'real time" PCR- See, 
for example, . WO/9746707A2, WO/974671 2A2, 
WO/974671 4A1, all published December 11, 1997, 
. incorporated herein by reference. 
[0003] Monitoring fluorescence each cycle of PCR 
initially involved the use of ethidium bromide. R. 
Higuchi, G. Dollinger, P. S. Walsh, and R. Griffith, 
"Simultaneous amplification and detection of specific 
DNA sequences," Bio/Technology 10:413-417. 1992; R. 
Higuchi, C. Fockler, G. . Dollinger, and R. Watson, 
"Kinetic PCR analysis: real time monitoring of DNA 
amplification reactions," Bio/Technology 11:1026-1030, 
1 993. In that system fluorescence is measured once per 
cycle as a relative measure of product concentration. 
Ethidium bromide detects double stranded DNA; if tem- 
plate is present fluorescence intensity increases with 
temperature cycling. Furthermore, the cycle number 
where an increase in fluorescence is first detected 
increases inversely proportionally to the log of the initial 
template concentration. Other fluorescent systems have 
been developed that are capable of providing additional 
data concerning the nucleic , acid concentration and 
sequence. 

[0004] While PCR is an invaluable molecular biol- 
ogy tool, the practical implementation of real time PCR 
techniques has lagged behind the conceptual promise. 
Currently available instrumentation does, not actually 
analyze data during PCR; it simply acquires the data for 



later analysis. After PCR has been completed, multiple 
manual steps are necessary to analyze the acquired 
data, and human judgment is typically required to pro- 
vide the analysis result. What is needed is a system for 
5 automating data acquisition and analysis so that no 
user intervention is required for reporting the analytical 
results. Thus, when the temperature cycling in a 
polymerase chain reaction amplification is complete, 
the system software is automatically triggered and the 
10 results, for example, the presence or absence of a given 
pathogen, are immediately displayed on screen. Algo- 
rithms for detection, quantification, and genotyping are 
needed. Moreover, initiation of the analysis algorithm 
can be implemented prior to completion of temperature 
>5 cycling. Data processing can occur during amplification 
and concomitant analysis results can be used to modify 
temperature cycling and to acquire additional data dur- 
ing the latter stages of the amplification procedure to 
optimize amplification protocol and data quality. 
20 [0005] A major problem in automating PCR data 
analysis is identification of baseline fluorescence. Back- 
ground fluorescence varies from reaction to reaction. 
Moreover, baseline drift, wherein fluorescence 
• increases or decreases without relation to amplification 
25 of nucleic acids in the sample, is a common occurrence. 
Prior attempts to automate amplification data analysis 
involved setting the baseline fluorescence as that meas- 
ured at one or more predetermined early cycle num- 
bers. This technique accounts for the variation in 
30 background fluorescence, but it does not compensate 
for baseline drift. Without compensation for baseline 
drift, automated amplification data analysis can easily 
provide both false negative and false positive results. 
[0006] Thus, one aspect of the present invention is 
35 directed to a method of determining the presence of a 
nucleic acid in a sample by using a fluorescent entity 
capable of detecting the nucleic acid and amplifying the 
nucleic acid in the presence of the fluorescent entity. A 
baseline fluorescence region is determined by analyz- 
es ing the fluorescence measurements of a number of 
amplification cycles, and the fluorescence measure- 
ments during specific amplication cycles are com- 
pared to the baseline fluorescence region to determine 
the presence or absence of the nucleic acid. In a pre- 
45 ferred embodiment, the baseline fluorescence region is 
determined by calculating the slope of a fluorescence 
intensity verses amplification cycle plot at each of the 
amplification cycles and choosing the fluorescence 
measurement of the amplification cycle with the slope 
so having an absolute value closest to zero. Preferably, the 
baseline fluorescence region is generated as the 
square root of the mean square error. 
[0007] In another embodiment, the baseline fluo- 
rescence region is determined and the fluorescence 
55 values are compared thereto after each amplification 
cycle. Thus, the presence of the nucleic acid sequence 
can be determined rnore quickly in samples containing • 
higher copy number. Furthermore, the remaining cycles 
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may be used to acquire other information concerning 
the nucleic acid sample, such as initial copy number and 
allelic data. 

[0008] In an additional embodiment,. the process is 
automated so that a user can prepare a biological sam-. 
pie and simply place it in a thermal cycler having a sen- 
sor for reporting fluorescence values as a function of 
cycle number and a processor programmed with an 
algorithm capable of processing the values and report- 
ing a positive or negative result. 
[0009] Additional features of the present invention 
will become apparent to those skilled in the art upon 
consideration of the following detailed description of 
preferred embodiments exemplifying the best mode of 
carrying out the invention as presently perceived. 

Rriof n^ri ption of the Drawings 

[001 0] 

Figs. 1A-L show a comparison of three fluores- 
cence monitoring schemes, dsDNA dye, exonucle- 
ase probe, and hybridization probe, for PCR 
amplification, wherein each scheme is illustrated 
before amplification (Figs. 1A-C) and after amplifi- 
cation (Figs. 1D-F), and fluorescence values are 
shown once during each cycle of PCR (Figs. 1G-I) 
and continuously during PCR (Figs. 1J-L). 
Fig. 2 is a graph illustrating logistic growth. 
Figs. 3A-F show a comparison of various cycle- 
verses-fluorescence curve types. 
Fig. 4 illustrates a sliding window analysis for deter- 
mining the slope of the fluorescence-verses-cycle 
number graph at each cycle. 
' Figs. 5A-B show typical fluorescence verses ampli- 
fication cycle graphs for (Fig. 5A) a negative sample 
and (Fig. 5B) a positive sample. 
Figs. 6A : C also shows typical amplification graphs 
wherein Fig. 6A shows fluorescence verses amplifi- 
cation cycle, Fig. 6B is the first derivative of fluores- 
cence verses amplification cycle, and Fig. 6C is the 
second derivative of fluorescence verses amplifica- 
tion cycle. 

Figs. 7-1 1 show the results for various samples 
wherein open white circles represent the fluores- 
cence measurement at each cycle, open black cir- 
cles represent the first derivatives, closed black 
circles represent second derivatives, large black 
circles connected by lines represent the points con- 
tributing to the baseline calculation, and the hori- 
zontal lines illustrate the baseline region. Figs. 7 
and 8 illustrate positive results, while Figs. 9-1 1 
illustrate negative results. 

nptaited Descript ion of the Invention 

[001 1] In describing and claiming the invention, the 
following terminology will be used in accordance with 



the definitions set forth below. 

[0012] As used herein, "nucleic acid," "DNA," and 
similar terms also include nucleic acid analogs, i.e., 
analogs having other than a phosphodiester backbone. 
5 For example, the so-called "peptide nucleic acids," 
which are known in the art and have peptide bonds 
instead of phosphodiester bonds in the backbone, are 
considered within the scope of.the present invention. 
[0013] As used herein, "fluorescence resonance 
w energy transfer pair- or "FRET pair" refers to a pair of 
fluorophores comprising a donor fluorophore and 
acceptor fluorophore, wherein the donor fluorophore is 
capable of transferring resonance energy to the accep- 
tor fluorophore. In other words the emission spectrum of 
15 the donor fluorophore overlaps the absorption spectrum 
of the acceptor fluorophore. In preferred fluorescence 
resonance energy transfer pairs, the absorption spec- 
trum of the donor fluorophore does not substantially 
overlap the absorption spectrum of the acceptor fluoro- 
20 phore. 

[0014] As used herein, "FRET oligonucleotide pair 
refers to a pair of oligonucleotides, each labeled wit a 
member of a fluorescent resonance energy transfer 
pair, wherein hybridization to complementary target 
25 nucleic acid sequences brings the fluorescent entities 
into a fluorescence resonance energy transfer relation- 
ship. 

[0015] ■ The present invention is directed to a 
method of analyzing a sample for the presence of a 
30 nucleic acid wherein the sample is amplified, preferably 
using PCR, in the presence of a fluorescent probe capa- 
ble of detecting the presence of the nucleic acid sample. 
A baseline region is determined by comparing the fluo- 
rescence at various amplification cycles, and the fluo- 
35 rescence at each of various amplification cycles is 
compared to the baseline region to determine whether 
the fluorescence measurements fall outside of that 
baseline region. 

[0016] Many different probes have recently become 
40 available for monitoring PCR. Although not sequence 
specific, double stranded DNA (dsDNA) specific dyes 
can be used in any amplification without the need for 
probe synthesis. Such dyes include ethidium bromide 
and SYBR™ Green I. With dsDNA dyes, product specif- 
45 icity can be increased by analysis of melting curves or 
by acquiring fluorescence at a high temperature where 
nonspecific products have melted. K. M. Ririe, R. P. 
Rasmussen, and C. TWittwer, "Product differentiation 
by analysis of DNA melting curves during the polymer- 
50 ase chain reaction," Anal. Biochem. 245-154-160, 
1997; T. B. Morrison, J. J. Weis, and C. T. Wittwer, 
"Quantification of low copy transcripts by continuous 
SYBR Green I monitoring during amplification," Bio- 
Techniques 24:954-962, 1998. 
55 [0017] Oligonucleotide probes can also be cova- 
lently labeled with fluorescent molecules. Hairpin prim- 
ers (Sunrise™ primers), hairpin probes (Molecular 
Beacons™) and exonuclease probes (TaqMan™) are 



XXSID: <EP_ 



EP 1 059 523 A2 



dual-labeled oligonucleotides that can be monitored 
during PCR. These probes depend on fluorescence 
quenching of a fluorophore by a quencher on the same 
oligonucleotide. Fluorescence increases when hybridi- 
zation or exonuclease hydrolysis occurs. 
[0018] A preferred design employs two oligonucle- 
otides, each labeled with a fluorescent probe. Hybridi- 
zation' of these oligonucleotides to a target nucleic acid 
brings the two fluorescent probes close together to 
allow resonance energy transfer to occur. C. T. Wittwer, 
M. G. Herrmann, A. A. Moss and R. P. Rasmussen, 
"Continuous fluorescence monitoring of rapid cycle 
DNA amplification," BioTechniques 22:130-138, 1997. 
These hybridization probes require only a single fluo- 
rescent label per probe and are easier to design and 
synthesize than dual labeled probes. Acceptable fluoro- 
phore pairs for use as fluorescent resonance energy 
transfer pairs are well known to those skilled in the art 
and include, but are not limited to, fluorescein/rhodam- 
ine phycoerythrin/Cy7, fluorescein/Cy5, fluores- 
cein/Cy5.5. fluorescein/LC Red 640, and fluorescein/LC 
Red 705. 

[0019] SYBR™ Green I, exonuclease probe and 
hybridization probe designs are shown in Figs. 1A-L. 
For each design, schematics both before (Figs. 1A-C) 
and after (Figs. 1 D-F) amplification are shown, as well 
as cycle verses fluorescence amplification plots of pos- 
itive and negative controls (Figs. 1G-I), and temperature 
verses fluorescence plots from continuous monitoring 
(Figs. 1J-L). SYBR Green I fluorescence increases as 
more dsDNA is made (Fig. 1A). Because the dye is not 
sequence specific, a negative control also increases in 
fluorescence during later cycles as primer dimers are 
formed (Fig. 1G). In Fig. 1B, dual-labeled fluores- 
cein/rhodamine probes are cleaved during polymerase 
extension by 5'-exonuclease activity, separating the 
fluorophores and increasing the fluorescein emission. 
The signal generated is cumulative and the fluores- 
cence continues to increase even after the amount of 
product has reached a plateau (Fig. 1 H). Fig. 1 C shows 
use of a FRET oligonucleotide pair wherein two probes 
hybridize next to each other, one labeled 3' with fluores- 
cein and the other labeled 5' with Cy5. As product accu- 
mulates during PCR, fluorescence energy transfer to 
Cy5 increases (Fig. 11). The fluorescence of hybridiza- 
tion probes decreases at high cycle number because of 
probe/product competition. 

[0020] Standard instruments for PCR complete 30 
cycles in about two to four hours. A preferred system is 
a rapid thermal cycling device using capillary tubes and 
hot. air temperature control. See, for example, United 
States Patent No. 5,455,175. herein incorporated by ref- 
erence. Because of the low heat capacity of air and the 
thin walls and high surface area of capillary tubes, small 
volume samples could be cycled quickly. The total 
amplification time for 30 cycles is reduced to 1 5 minutes 
with excellent results. 

[0021 ] The use of capillaries with forced air heating 



allows precise control of sample temperature at a speed 
not possible with other designs. For example, sample 
temperature verses time plots in capillaries show sharp 
spikes at denatu ration and annealing temperatures, 
; whereas several seconds are required for all of the sam- 
ple to reach equilibrium in conical plastic tubes. C. T. 
Wittwer G. B. Reed, and K. M. Ririe, "Rapid cycle DNA 
amplification," in K. Mullis, F. Ferre, and R. Gibbs (Eds.), 
"The polymerase chain reaction," Springer-Verlag. 
J0 Deerfield Beach. FL. pp. 174-181, 1994; C. T. Wittwer, 
B C. Marshall, G. B. Reed, and J. L. Cherry, "Rapid 
cycle allele-specific amplification: studies with the cystic 
fibrosis delta F508 locus," Clin. Chem., 39:804-809, 
1993. Rapid temperature cycling with minimal anneal- 
, 5 ing and denaturation times improves quantitative PCR 
and increases the discrimination of allele specific ampli- 
fication. J. H. Weis, S. S. Tan, B. K. Martin, and C. T. 
Wittwer, "Detection of rare mRNA species via quantita- 
tive RT-PCR, Trends in Genetics," 8:263-4, 1992; S. S. 
20 Tan and J. H. Weis, "Development of a sensitive reverse 
transcriptase PCR assay, RT-RPCR, utilizing rapid cycle 
times," PCR Meth. and Appl. 2:137-143, 1992. Rapid 
cycling for cycle sequencing reduces sequencing arti- 
facts and minimizes "shadow banding" in dinucleotide 
25 repeat amplifications. H. Swerdlow, K. Dew-Jager, and 
R F Gesteland, "Rapid cycle sequencing in an air ther- 
mal cycler," BioTechniques 15:512-519. 1993; S. J. 
Odelberg and R. White, "A method for accurate amplifi- 
cation of polymorphic CA-repeat sequences." PCR 
30 Meth. Appl. 3:7-12. 1993. For long PCR, yield is 
improved when the sample is exposed as little as possi- 
ble to high denaturation temperatures. C. E. Gustafson, 
R. A. Aim and T. J. Trust, "Effect of heat denaturation of 
target DNA on the PCR amplification," Gene 23:241- 
3 5 244 1993. The RapidCycler™, developed by Idaho 
Technology, is an example of a rapid thermal cycling 
device. The LightCycler™ is a rapid temperature cycler 
with a fluorimeter, wherein light emitting diodes are 
used for excitation and photodiodes are used for detec- 
40 tion. 

[0022] The present invention is directed to methods 
for automating detection nucleic acids with real time 
PCR. While these algorithms may be applied to any 
. amplification system, it is preferred to integrate these 
45 algorithms into the LightCycler™ platform. These analy- 
sis routines are triggered by the completion of rapid 
thermal cycling for "hands off amplification, analysis, 
and final results presentation in a total of less than 15 
min. The analysis routines take from <1 second for 
so detection and quantification to <10 seconds for geno- 
typing LabView (National Instruments, Austin, TX), a 
graphical programming language, is preferred for Light- 
Cycler™ instrument control. The LightCycler™ is a PC- 
based instrument. 
55 [0023] Perhaps the most basic analysis of real time 
PCR data is a judgement of whether a targeted nucleic 
acid is present. If the nucleic acid is present, further 
quantification and genotyping may take place. In many 
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cases, a yes/no judgement is all that is needed. For 
example, one may want to determine whether E. coli 
0157:H7 is in a sample of hamburger, whether anthrax 
is on a swab from a soldier; or whether hepatitis C is in 
a unit of blood. Real time PCR can improve yes/no 
detection over end point PCR assays because fluores- 
cence is acquired at each cycle. 

[0024] Inspection of cycle verses fluorescence data 
from positive and negative real time PCR runs (see Fig. 
1H and 11) suggests that discrimination is simple. The 
positive samples increase with cycle number while the 
negative samples remain at baseline. A trained 
observer expects positive samples to follow an S-shape 
curve, beginning with a baseline, followed by an expo- 
nential segment, and finishing with a plateau. The 
expected curve is similar to the logistic model for popu- 
lation growth, where the rate of growth is proportional to 
both the population size y and to the difference L-y, 
where L is the maximum population that can be sup- 
ported. For small y, growth is exponential; but as y nears 
L. the growth rate approaches zero. An example of 
logistic growth is shown in Fig. 2. 
[0025] Although intuitively simple, accurately dis- 
criminating between positive and negative samples is 
not easy in practice. The simplest approach is to set a 
horizontal fluorescence threshold as a discriminator 
between positive and negative samples. This works 
best with a stable baseline (between and within sam- 
ples) and a known fluorescence intensity that correlates 
with "positive." Although this method will work on obvi- 
ous samples (e.g., Fig. 1 H and 11), a more robust algo- 
rithm is desired that will work under a wider variety of 
conditions. For example, the baseline may drift and the 
fluorescence intensity may vary greatly between differ- 
ent samples and probe techniques. Thus, the present 
invention is directed to a method that will: (1) automati- 
cally identify the baseline, (2) use the baseline variance 
to establish a confidence region, and (3) call each sam- 
ple positive or negative based on the relationship of the 
confidence region to the fluorescence data. 
[0026] Figs. 3A-F display various types of amplifica- 
tion curves, all of which have been observed in LightCy- 
cler™ runs. Figs. 3A and B show curves from samples 
that are negative with no template present. The fluores- 
cence scales in Figs. 3A and B are magnified (com- 
pared to Figs. 3C-F) to demonstrate the baseline drift 
and to provide algorithms capable of being independent 
of the fluorescence intensity. There, is always some 
baseline drift during cycling. This drift usually is greatest 
at the beginning of cycling but later levels off. and may 
be either downward (Fig. 3A) or upward (Fig. 3B). This 
baseline drift of negative reactions must be distin- 
guished from positive reactions of either low copy num- 
bers (Fig. 3C) or high copy numbers (Fig. 3D) of starting 
template. The method needs to work with various probe 
designs, including exonuclease (Fig. 3E) and hybridiza- 
tion (Fig. 3F) probes. 

[0027] Automatic identification of the background is 



surprisingly difficult. In prior art methods, the baseline is 
determined as a function of measured fluorescence at a 
fixed range of cycles near the beginning of amplifica- 
tion However, selection of a fixed range of cycles is not 
5 adequate because both downward drift (Fig. 3A) and 
high copy (Fig. 3D) amplifications may be incorrectly 
called. 

[0028] In the present invention, the background is 
identified by analyzing the fluorescent measurements 
w over a wide range of amplification cycles. Preferably, the 
background is identified by selecting the sliding window 
(Fig. 4) with the shallowest slope. That is, calculate the 
slope at each cycle by liner regression of the local 
neighborhood (for example, a 7 point sliding window). 
is The window with the slope of lowest absolute value 
(least difference from zero) defines the background, 
region. Once the background region has been identi- 
fied, the variation of these background points about 
their regression line (the square root of the mean 
20 square error) is multiplied by a constant to determine a 
confidence band. This confidence band will have a 
slope near zero and is extrapolated across all cycles. If 
the fluorescence of the last cycle is within the confi- 
dence band it is negative, if it is outside the band it is 
25 positive. Figs. 5A-B demonstrate both cases. . 

[0029] This algorithm should work well in most 
cases. However, with the high copy fluorescence curve 
type (Fig. 3D), the shallowest slope might be found at 
early cycles (resulting in a correct positive call) or at late 
30 cycles (resulting in an incorrect negative call). This 
exception may be handled by analyzing the curve 
shape. In a well-behaved amplification, the expected 
amplification curve shape is ordered by cycle number 
as follows: 



1. Minimum fluorescence 

2. Maximum second derivative (F") 

3. Maximum first derivative (F') 

4. Minimum second derivative (F") 
40 5. Maximum fluorescence 

This gives the characteristic S-curve shape we expect 
during PCR (Fig. 6A). The maximum slope (first deriva- 
tive) is obtained from the sliding window analysis 
45 already performed for background identification. Prefer- 
ably, the second derivatives are calculated by a 3-point 
sliding window liner regression of the first derivatives. If 
the curve shape is well behaved (that is, if looking at a 
graph of Figs. 6A-C, and reading from lowest to highest 
so cycle number, the features occur in the order listed 
above), ten the background is only selected from sliding 
windows centered at cycle numbers less tan the second 
derivative maximum. This solves the potential analysis 
problem wit Fig. 3D. In other preferred embodiments, 
55 cycle numbers less than the first derivative maximum or 
cycle numbers less than the second derivative minimum 
may be used. It will be further understood that any cycle 
number between the second derivative maximum and 
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the second derivative minimum is a suitable cutoff cycle 
for use with this technique and is within the scope of this 
invention. 

[0030] Another method is to compare the cycle with 
the greatest "fluorescence (which is not necessarily the 
last cycle) to the confidence band. This is especially 
suited for hybridization probes tat may decrease in fluo- 
rescence with extensive cycling, such as seen in Fig. 
3F. The cycle with the greatest fluorescence only should 
be used when the curve shape is well behaved, in order 
to prevent false positive calls with downward drifts, such 
as shown in Fig. 3A. 

[0031] The variables to optimize for automatic 
detection are: 1) the window size for the first derivative 
estimate, 2) the window size for the second derivative 
estimate, and 3) the confidence band factor. A reasona- 
ble value for the first derivative window size is 7, 
although 3, 5, 9, and 1 1 are also quite useful. For the 
second derivative the preferred window size is 3, but 5, 
and 7 have also proven to be useful values. A preferred 
confidence band factor is 20. As the first derivative win- 
dow size increases the variance estimate is more accu- 
rate, but the edge cycles (beginning and ending) are 
lost. 

[0032] This algorithm is best understood by refer- 
ring to the fluorescence verses cycle test result plot 
shown in Figs. 7-1 1 . The input data consist of one fluo- 
rescence value for each cycle of amplification, shown as 
the closed white circles. Let this equal array Yi, where i 
is the cycle number and N is the total number of cycles. 
The detection criteria are: 

A = the number of fluorescence values used to 
determine the first derivatives. 

It is convenient to use odd numbers, so that the 
first derivatives correspond to integer cycle num- 
bers. As discussed above, reasonable values 
include 3, 5, 7, 9, and i 1. Preferably, 7 is used as 
■ the first derivative window size. 

B = the number of first derivative values used to 
determine the second derivatives. 

Again, it is convenient to use odd numbers, so 
that the second derivative values also correspond 
to integer cycle numbers. Reasonable values 
include 3, 5, and 7, with 3 being the preferred value. 
C = the confidence band factor. 

This factor determines the confidence band by 
multiplying it by a variance measure, preferably the 
square root of the mean square error. 

[0033] The first step is to calculate the first and sec- 
ond derivatives. Although there are many ways to 
accomplish this, a preferred method is to determine the 
first derivatives as the slope of a linear regression line 
through A points, and assigning the value to the central 
cycle number. Some cycles on either edge cannot be 
assigned first derivatives, but first derivatives can be 
provided for cycles (A + 1)/2 through N- (A - 1)/2 . Sim- 



ilarly, the second derivatives are calculated as the slope 
of the first derivative points and assigned to cycles 
(A + 1 )/2 + (B - 1)/2 through [N - (A - 1)/2] - (B - 1)/2 . 
Calculation of the first anc second derivatives provide 
5 arrays Y'i and Y"i, with some edge values missing. In 
Fig. 7, the first and second derivatives are displayed as 
open black circles and closed black circles, respectively. 
[0034] The next step is to determine whether the 
fluorescence curve has a well-behaved shape. As dis- 
io cussed above, the well-behaved shape occurs when the 
cycles with minimum fluorescence, maximum second 
derivative, maximum first derivative, minimum second 
derivative, and maximum fluorescence occur in that 
order, from law to high cycle number. 
15 [0035] The baseline is then determined. If the fluo- 
rescence curve does not have the expected shape, the 
cycle whose first derivative is closest to zero is used. If 
the fluorescence curve has a well-behaved shape, the 
cycle whose first derivative is closest to zero chosen 
20 from among all cycles prior to the cycle with the maxi- 
mum second derivative (again, any cycle between the 
maximum second derivative and the minimum second 
derivative may also be used as the cutoff cycle number). 
The baseline is drawn through the fluorescence value of 
25 the chosen cycle with a slope of its first derivative. In 
Fig. 7, the A points contributing to the first derivative cal- 
culation for the baseline are displayed as large black 
dots connected by a line. 

[0036] The next step is to determine the test point 
30 cycle, that is, the cycle used to compare against the 
baseline for determining a positive or negative result. It 
the curve is not well-behaved, the test point is the last 
cycle. If the fluorescence curve is well-behaved, the test 
point is the cycle with fluorescence farthest from the 
35 baseline. The test point fluorescence of a negative sam- 
ple can be predicted as the intersection of the baseline 
with the test point cycle. 

. [0037] Next, a confidence interval can be deter- 
mined about the predicted negative test point. Prefera- 
40 bly, this is done by finding the square root of the mean 
square error about the baseline of A points used to 
determine the 1 baseline. This is multiplied by C. The 
product is added to the predicted negative test point to 
get the upper fluorescence limit of the confidence inter- 
45 val and is subtracted from the predicted negative test 
point to get the lower limit of the confidence interval. 
These limits are shown on Fig. 7 as two solid horizontal 
lines. 

[0038] The final step is to declare the sample posi- 
50 five or negative. If the test point fluorescence is outside 
of the confidence interval, the sample is positive. If it is 
within the interval, the sample is negative. Figs. 7 and 8 
are samples which are positive, while Figs. 9-11 are 
negative samples. 
55 [0039] Although the invention has been described 
in detail with reference to preferred embodiments, vari- 
ations and modifications exist within the scope and spirit 
of the invention as described and defined in the follow- 
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Claims 

1 . A method for determining the presence of a nucleic s 
acid comprising the steps of 

providing a fluorescent entity capable of indi- 
cating the presence of the nucleic acid and 
capable of providing a signal related to the 10 
quantity of the nucleic acid, 
amplifying the nucleic acid through an initial 
plurality of amplification cycles in the presence 
of the fluorescent entity, 

measuring the fluorescence of the fluorescent 15 
entity during each of the initial plurality of ampli- 
fication cycles, 

analyzing the measured fluorescence to iden- 
tify amplification cycles for use in establishing a 
baseline fluorescence region, and 20 
ascertaining whether the fluorescence meas- 
urement during any of the initial plurality of 
amplification cycles is outside the baseline flu- 
orescence region. 

2i 

2. The method of claim 1 wherein the cycles used in 
establishing the baseline fluorescence region are 
identified by calculating slopes of a fluorescence 
verses amplification cycle plot at the plurality of 
amplification cycles. 

3. The method of claim 2 wherein the slopes are cal- 
culated by linear regression of each local neighbor- 
hood. 

s 

4. The method of claim 2 wherein the cycles used in 
establishing the baseline fluorescence region com- 
prise an interval of cycles that includes the amplifi- 
cation cycle with the slope having an absolute value 
closest to zero 

5. The method of claim 4 wherein the amplification 
cycle with the slope having the absolute value clos- 
est to zero comprises a central cycle of the- interval 
of cycles used to establish the baseline fluores- 
cence region. 

6. The method of claim 4 wherein the baseline fluores- 
cence is centered on a line passing through the 
measured fluorescence of the amplification cycle 
with the slope having the absolute value closest to 
zero, the line further comprising a slope equal to the 
slope having the absolute value closest to zero, and 
the baseline fluorescence region including an area 
. on at least one side of the line defined by a variance 
measure of the fluorescence values. 

7. The method of claim 6 wherein the variance meas- 



ure is calculated from the initial amplification cycles 
used in establishing the baseline fluorescence. 

The method of claim 1 wherein the fluorescent 
entity comprises a FRET oligonucleotide pair. 

The method of claim 1 further comprising the steps 



generating a well-behaved amplification curve, 
determining a cutoff cycle, 
calculating the slope of the fluorescence 
verses cycle plot for each of the initial amplifi- 
cation cycles prior to the cutoff cycle, and 
choosing the fluorescence measurement of the 
amplification cycle with the slope having an 
absolute value closest to zero. 

10. The method of claim 9 wherein the cutoff cycle is 
selected from the group consisting of a maximum 
second derivative, a maximum first derivative, and a 
minimum second derivative of a fluorescence 
verses cycle curve. 

11. The method of claim 1 further comprising the steps 



generating a well-behaved amplification curve, 
determining the cycle having the maximum flu- 
orescence, and 

determining whether the maximum fluores- 
cence is outside the baseline fluorescence 
region. 

12. The method of claim 1 further comprising the steps 
of 

generating an amplification curve which is not 
well behaved, and 

determining if the fluorescence of a last cycle 
tested is outside the baseline fluorescence 
region. 

13. The method of claim 1 further comprising the steps 



assessing whether a positive result is indicated 
by determining whether the fluorescence 
measurement during any of the of amplification 
cycles is outside the baseline fluorescence 
region, and 

continuing amplification through a plurality of 
additional amplification cycles and repeating 
the measuring, determining, ascertaining, and 
assessing steps after each additional amplifi- 
cation cycle until either the positive result is 
obtained or a maximum cycle number is 
reached. 
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14. The method of claim 13 further comprising the 
steps of 

performing a second plurality additional ampli- 
fication cycles after the positive result is 5 
obtained, and 

analyzing the nucleic acid for additional infor- 
mation. 

15. The method of claim 13 wherein the nucleic acid is to 
further analyzed to determine the presence of a 
particular allele. 

16. The method of claim 13 wherein the measured flu- 
orescence is further analyzed to determine initial is 
concentration. 

17. An automated method for determining the presence 
of a nucleic acid comprising the steps of 

20 

placing a sample into a container containing a 
fluorescent entity capable of indicating the 
presence of the nucleic acid and capable of 
providing a signal related to the quantity of the 
nucleic acid, 25 
placing the container into a device for amplify- 
ing the nucleic acid through a plurality of ampli- 
fication cycles in the presence of the 
fluorescent entity, 

measuring the fluorescence of the fluorescent 30 
entity during each of the plurality of amplifica- 
tion cycles, 

determining a baseline fluorescence region by 
analyzing the fluorescence at various amplifi- 
. cation cycles, and 35 
outputting a positive result if the fluorescence 
measurement during any of the plurality of 
amplification cycles is outside the baseline flu- 
orescence region. 



50 
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(54) Automated real-time analysis of nucleic acid amplification 



(57) A method is described for analyzing a sample 
for the presence of a nucleic acid wherein the sample is 
amplified, preferably using PCR,. in the presence of a 
fluorescent probe capable of detecting the presence of 
the nucleic acid sample. A baseline region is determined 
by comparing the fluorescence at various amplification 



cycles, and the fluorescence at each of various amplifi- 
cation cycles is compared to the baseline region to de- 
termine whether the fluorescence measurements fall 
outside of that baseline region. 
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